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Abstract 



Double integrals that represent matrix elements of the power and logarithmic potentials 
(resp. \x — y\ a and \n\x — y\ ) in the Legendre polynomial basis are found in a closed 
form. Several proofs are given, which involve different special functions and identities. In 
particular, a connection of the new formulae and Whipple's hypergeometric summation 
formula is shown. 
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§ 1. Introduction 

In this paper we evaluate the integrals 

1 1 

B a (m,n) = J J \x - y\ a P m (x) P n (y) dx dy (Re a > -1) (1.1) 

-i-i 

(matrix elements of the Bessel potentials) and 

i 

B a (m,n) = J J \n\x-y\P m (x)P n (y)dxdy. (1.2) 



B'(m, n) = — — 



da 



o=0 



-i-i 



Here P n are the Legendre polynomials. Notation and relevant facts regarding orthogonal 
polynomials and hypergeometric series are collected in Appendix. 

Sometimes it is more convenient to deal with the shifted Legendre polynomials (jA.6|) . 
Define 

i i 

B a (m,n) = J J \x - y\ a pm(x) p n (y) dx dy (Rea > -1) (1.3) 
o o 

and 

i i 

B'(m,n) = J J In \x - y\ p m {x) p n (y) dx dy. (1.4) 
o o 



It is readily seen that 



B a (m,n) = 2 a+2 B a (m,n), 

B'(m,n) = 4 ( B'(m,n) + 5 m0 5 n0 In 2 



Because of the symmetry of the integrands — see Appendix A. 1.3. — all the above integrals 
vanish when m — n is odd. In fact, we obtain more detailed results, which are non-trivial in 
the case of odd m — n as well, — evaluation of the 'halves' of the above integrals. Instead 
of the interals over the square 

-1 < x < 1, -1 < y < 1 ( resp. < x < 1, < y < 1 ), 
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consider the integrals of the same integrands over the triangle 

< y < x < 1 (resp. < y < x < 1), 
and denote them by L a (m,n), L'(m,n), L a (m,n), L'(m,n) respectively. Then 

L(m,n) = (-ir- n L(n,m); 

(1.5) 

B(m,n) = 2L(m,n) if m — n is even, 

with L being any of L a , L', L a , U and the corresponding B. Yet we formulate theorems 
on evaluation of the S's separately, because there are compact proofs that do not involve 
evaluation of the L's. 

In view of the symmetry (|1.5|) . it is enough to bring the results assuming m > n. Denote 
d = m — n. Set 

R a (m,n) = _=L J id Z a ^t { 7~n, d (1-6) 

and 

R'(m,n) = — r— r. (1.7) 

V ' ; (m + n){m + n + 2)(d? - 1) V ; 

Notice that if d > then the fraction R a (m,n) can be reduced by the common factor 
(a + 1). 

Theorem 1. If in — n > is even, then 

B a (m,n) = 2R a (m,n). (1.8) 

Theorem 2. (i) If m — n is even and (m,n) ^ (0,0), then 

B'(m,n) = 2R'(m,n). (1.9) 

(ii) Special case m = n = : 

B'(0,0) = -3/2. (1.10) 

Correspondingly, for even nonnegative d = m — n, 

B a (m,n) = 2 a+3 R a (m,n), (1.11) 

In particular, 

2^+3 

B a (0,0) 



(a+l)(a + 2)" 
Also, assuming m = n even and (m, n) 7^ (0, 0), we have 

B'(m,n) = 8R'(m,n). (1.12) 

The special case is 

B'(0, 0) = 4 In 2 - 6 = -3.2274 .... (1.13) 



Theorem 3. If m — n > then 

L a (m,n) = R a (m,n). 

Theorem 4. (i) // m — n > 2 or m = n ^ 0, i/ien 

L'(m,n) = R'(m,n). 

(ii) Special case m — n = 1 : 



( rxodd _ _|_ _ _J_ 



where 



L'(n + l,n) = 
g = (2n+l)(2n + 3) = (to + n)(m + n + 2), 



;i.i4) 



;i.i5) 



;i.i6) 



and -^2n+i ^ s ^ e n + °dd harmonic sum, 



TT 



odd 
2n+l 



1 1 
1 + - + ... + -. 

3 2n + 1 



The following numerical table supports Theorem 4 and includes the omitted case m 
Index m enumerates the rows and n the columns. 



(1.17) 

71 = 0. 



L'(m, n) 



m,n=0...3 





- -3 


5 


1 


1 




~T 


36 


24 


120 




-5 


-1 


61 


1 




36 


"8" 


900 


72 




1 


-61 


-1 


527 




24 


900 


24 


14700 




-1 


1 


-527 


1 




- 120 


72 


14700 


120 



The paper is organized as follows. 

Section 2 contains a one-page proof of Theorem 2, due to M. Rahman jllj . 

Section 3 contains a short proof of Theorem 1 proposed by R. Askey. It involves three 
interesting formulae from the theory of special functions. The proof is simple because nasty 
details are worked out in the cited results. 

In Section 4 we prove Theorem 3 by transforming the integral L a (m,n) to a certain 
terminating hypergeometric series, which can be evaluated by a formula due to Whipple. 
Moreover, through this connection we prove the terminating case of Whipple's formula. 
(Essense: our integral is a rational function with known zeros and poles.) A proof of 
Whipple's formula in the general case is then obtained using a standard function-theoretic 
argument. 

Finally, in Section 5 we give an independent, self-contained, and elementary proof of 
Theorem 4, which had been largely computer-aided. 
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This work stems from my study of the integral B'(m, n) that was originally motivated by 
numerical solution of integral equations on polygons. Formula (jl.l2|) was published without 
proof in |13j . The following problem about a generalization of the integral B'(m,n) is 
important in those applications. 

Problem. Propose an efficient numerical procedure for evaluation of the integral 

1 1 

B'(m, n; A) = J J \n(Xx - y) p m (x) p n (y) dx dy, 
o o 

where X is a complex parameter. The method must be numerically stable and suitable for 
large values of m and n, and use the ordinary floating point arithmetics. 

The B'(m, n; A) are rational functions of A, if m — n > 2. Their evaluation using computer 
algebra with exact rational arithmetics is not what we wish, since as |m — n\ grows, the 
expressions become huge and the calculation slow. 

We conclude this section by deriving Theorem 2 from Theorem 1. 

Case n = d = stands alone, but it is trivial. Cases (i) d ^ and (ii) d = , n > 
need separate treatment, but in both cases B a is a rational function of a with a simple 
zero at a = 0. It remains to factor out the a and to evaluate the remaining quotient at 
a = 0. In case (i) the a is contained in (—a — 1)^. We have: as a — > 

B a (m,n) ^ 2(j) w (l) d _ 2 2j (| + 1) 

a ~* (f + 2) n (2) d+1 (f +n + l)(f + n)(d - l)d(d+ l)(d + 2)' 

Simplifying, we obtain ([1.9j) . In case (ii) the factor a comes from (d — a/2) n = (— a/2) n . 
When a — > 0, 

a- 1 S a (n,n) -> -(l)„_i/(2 • (2) n ) = -l(2n(n + l))" 1 , 
in agreement with ()1.9|) . 
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§ 2. M. Rahman's proof of Theorem 2 

In this proof, the double integral is turned to a sum of two ordinary integrals of a Legendre 
polynomial and a Legendre function of the 2nd kind, for which the values are known. We 
work with polynomials P n , rather than with the shifted Legendre polynomials, so the target 
is formula (|1.12j) . 

Lemma 2.1. Denote 

2n + 1 r 1 

F n (x) = — ^— jf_ i P n (y) In \x - y\ dy. (2.1) 

Then 

Fn (x) = l - p.v. / r» + i(v) - Pn-i(v) dy _ {2 2) 

2 J x — y 

Proof. Replacing In \x — y\ in (|2.1|) by ln(a; — y) doesnt't affect real part of the integral, 
whatever branch of the logarithm is taken. Move the integration path to the upper complex 
half-plane to obtain a regular integrand, and integrate by parts, using (jA.5|) . The boundary 
terms vanish, since [P n +i — -Pn-i](±l) = 0. Putting the integration path back on the real 
line, we obtain fl2.2(l . □ 

The integral in ()2.2|) is the difference of two Legendre functions of the second kind, Q n -i(x) 
and Q n {x) 0,3.6(29). Therefore, 



2n + 1 
2 

where 



B'(m, n) = J 1 P m {x) F n {x) dx = I(m, n - 1) - I(m, n + 1), (2.3) 



fi 

I(m,k) = I P m {x)Q k (x)dx. 



-i 

The integrals I(m,k) are evaluated in [U], 3.12 (17). The result is: 

(i) if m — k is even, then I(m, k) = ; 

(ii) if m — k is odd, then 



I(m,k) = — r. (2.4) 

v ' ; (m- k)(k + m + l) K J 



Assuming that m — n = d is even, substituion of ()2.4|) to ()2.3j) gives 
2n + 1 



B'(m, n) 



2 v ' ' (d+l)(m + n) (d- 1) (m + n + 2)' 
and p.!2|l easily follows. 
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§ 3. R. Askey's proof of Theorem 1 

In response to my question regarding formula posted in the online forum on special 

functions, R. Askey |3] wrote: 1 

There is a long paper by Polya and Szego in the early 1930s, see either Polya' s collected 
papers or those of Szego, in which they find an expansion of \x — y\ a as a sum of products 
of ultraspherical polynomials. Then the integral 11 . J)) can be evaluated as a single hyperge- 
ometric series by using the formula for the integral of an ultraspherical polynomial times a 
Legendre polynomial. This is a special case of a formula of Gegenbauer which I have used 
frequently. There are a number of relatively simple derivations of this formula. One of the 
easiest is the one Gasper and I gave in the paper we wrote for the meeting celebrating de 
Branges's proof of the Bieberbach conjecture. It just uses the generating function of ultra- 
spherical polynomials, differentiation, and orthogonality to get the result. Then it is just 
a matter of seeing if the hyper geometric series which comes from these calculations can be 
summed. Sadov claimed a specific formula is true, so the series can be summed. 

This section supplies details of Askey's approach. 

Step 1. Decomposition. Ultraspherical (or Gegenbauer's) polynomials p( u '(x) are 
defined by the generating function 



(1 - 2xt +t 2 )-" = Y. p n\x)t n - (3.1) 

n=0 

The decomposition of \x — y\ a due to Polya and Szego is jTUl S. 27] 

OO / • v 

r = M 1 £(l+ l) pM(x)Pjr\y), (3.2) 



\x 



3=0 



where v = —a/2 and 



rU-u) r(i + v ) 

M x = Mi(a) = V2 -. (3.3) 



In f!3.2j) . a > 0. Convergence in the r.h.s. follows from the estimation 

Pt"\x) = Oif- 1 ), -1 < x < 1. 
Step 2. Integration. The formula due to Gegenbauer reads |T| 



P^(x)P^ 2s (x)dx = frk-'fr 1/2)s . (3.4) 

Sl { L / Z )j-s+l 



Assume m — n = d > and even. Integral of the j-th term in (j3.2j) times P m (x) P n {y) may 
be nonzero only when j = m + 2k, k > 0. We have 



(prv), = (v) ^ 2 - 1/2)t 

(fffe), P„(y)) 



m+k+l 
+k+d/2 \y — l/2)fe+<i/2 



(fc + d/2)! (1/2) 



n+k+d/2+l 



1 Citation is slightly abridged. 
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Write the obtained sum in the standard hypergeometric form (jB.5|) . with a prefactor, using 
the identity 

i _ 2j_ = / + m\ ((m + v)/2 + l) k 



a V v) ((m + u)/2) k 
and the identities of the type ()B.2|) . for example, 

(v)n+k+d/2 = {v) n+d/2 {v + n + d/2) k . 

Also recall that n + d/2 = m — d/2. The result of integration takes the form 



B a (m,n) = Ml M 2 - 5 F 4 ( \ ' 2 ' 2 (3.5) 



with 

M = (1 + m/t/) (i/) m (i/) m -d/ 2 (y ~ l/2)d/ 2 , , 

(l/2) m+1 (l) d/2 (1/2) 

m— rf/2+1 

Step 3. Summation. The above series 5F4 can be summed by the limiting case (jB.llJ) of 
Dougall's formula for a well-poised terminating jFq. The result is a ratio of products of T 
functions. Arguments of four T's in the numerator are 

m + 3/2, rf/2 + 1, m+(3-d)/2, a + 2. (3.7) 

In the denominator we have four T's with arguments 

m + 1- a/2, (3 + a)/2, m + 2 + (a-d)/2, (3 + a + d)/2. (3.8) 

In total, we obtain an expression for B a (m, n) in the form of a fraction with 10 factors in the 
numerator and 8 factors in the denominator. To simplify, we break down the triple product 
()3.5|) into 5 groups and write 



BJm,n) = III n 2 IT3 n 4 IT 



5- 



In the expressions for the ITs below, we show which little factors come from which big ones. 
The notation [j, k] stands for the fc-th term in the numerator or denominator of the j-th big 
factor: j = 1 refers to M\, j = 2 to M 2 ; and j = 3 to 5F4. Denominator terms are shown 
in bold. For example, [2,1] is the first term in the numerator of ([3.6)1 . that is, (1 + mjv). 
One cancelation [3, 2]/[2, 2] = 1 occurs immediately and those terms are not attributed to 
any group. The grouping of other terms and the results of simplification are as follows 



ri! = 


[1,2] [2,1] [2,2]/[3,l] = 1, 






n 2 = 


[3,1]/[1,1] [2,1] = T(m + 3/2), 






n 3 = 


[l,l]/[3,2] = 2/(a + l), 






n 4 = 


[3, 4]/ [3, 3] [3,4] = 2 2+a+d n- 1/2 


\a + 2) d+1 (2 + 2±4) n 


-1 


n 5 = 


[2,2] [2,3] [3,3] (-a/2) m _ d/2 ( 


-(a + l)/2) d/2 r(m+^) 


[2,3] 


,l/2)m-d/2+l 





The part r(. . .) / (l/2) m _ d / 2+ i in II5 cancels with n 1 I 2 from II4. It remains to apply the 
duplication formula (|B.4j) with a = —a/2 , n = d/2, and (jl.ll)) follows. ■ 
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Limiting case: another proof of Theorem 2 



A proof of Theorem 2, different from those given in Sections 2 and 5, can be obtained as 
a special limiting case of the above argument. It follows the same lines, but instead of 
the ultraspherical polynomials we encounter Chebyshev's polynomials T n , where analogues 
of (|3.2|) and ()3.4|) are rather elementary. Chebyshev's polynomials show up here, because 
(2/n) T n (x) (n > 0) and the generating function for T n 's is 



d v PP{x 



u=0 



In (1 - 2 x t + t 2 ) = -2 >T n- 1 T n (x) f 



71=1 



The role of (|3.2j) is played by the formula 

oo 

ln|x-y| = -ki2 - 2j2r 1 T j {x)T j {y), 
which can be identified with the well known Fourier expansion 

oo 

-m|2sin(0/2)| = £ 1 ™s{j6) 
i=i 

using the definition of Chebyshev's polynomials Tj(cos8) = cos(j6). 
The analogue of ()3.4|) reads (for j > 0) 

2 x ^ , x , (« + l)j-2 a -l 



Pj _ 2s (x) dx 



J _{ Jy ' J " ' 2(s-l/2) i _ 2a+2 - 

Step 3 again amounts to evaluation of the particular case of the series ()3.6|) with v = 
using Dougall's formula. (It seems attractive to find its elementary evaluation in this case.) 

§ 4. Proof of Theorem 3 via Whipple's formula 

4.1. Reduction of integral L a (m,n) to a terminating 3F2 series 

Set 



A a (n, x) = (x -y) a p n (y)dy, (4.1) 
J 

then 

L a (m,n) = j p m (x) A a (n, x) dx. (4.2) 







Lemma 4.1. 

x a+1 f—n, n + 1 



Aafaz) = -— 2 F X ' ;x . (4.3) 
a + 1 V a + 2 / 

In particular, 

A a (n, 1) = { 7^ ( ~ a)n . (4-4) 

(a + l)n+l 
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Proof. Apply the Rodrigues formula (jA.7|) and integrate by parts n times. Assume a > n, 
then the boundary terms vanish and we get 

} (- a ) (-l)n 

A a (n, x) = j 1 ^ ) (x - y) a ~ n y n (1 - y) n dy. 
o 

Setting y = tx leads to the integral (jB.6|) for 2 i*i with parameters a = —n, b = n + 1, 
c = a + 2, and the prefactor 



a+l 



(-l) n (-q) n r(w + l)r(q-w + 1) 
n! T(a + 2) 



which reduces to x a+1 /(a + 1) . The assumption a > n can be dropped because both parts 
of (|4.3|) are analytic in a. □ 



Lemma 4.2. 



ft \ (-a-l) m „ ( ~n, n+1, a + 2 \ 

L a {m,n) = - — — 3^2 . i q , o • 45 

(a + l) m+2 \ a + m + 3, a + 2 - m J 

Proof. Substituting (JOI) in (jOjl . we find 

L a {m,n) = > ,, J , i / x + +J p m (x)dx. 

Since p m {x) = (— l) m p m (l — x), the integral in the common term is 

(-ir /(m, a + 1 + j, i) = , ( ~°:|~ j)m , 

(a + 2 + j) m+ i 

according to (|4.4|) . Using the identities (cf. (|B.2j) ) 

(a + l) i+ i (a + j + 2) m+1 = (a + l) m+2 (a + m + 3)j 

and 

{—01 — J ~ l) m = {-O! - l)r 



(a — m + 2)j 
the common term is transformed into 

(-a-l) m (-n) j (n + l) j (a + 2) j 



(a + l) m+2 (a + m + 3)j (a + 2 - m)j ' 
so dHU) follows. □ 

Remark. The symmetry (|1.5jl apparent in the l.h.s. of (|4.5jl is hidden in the r.h.s., where 
it is a particular case of Thomae's formula (|B.8|) . 
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4.2. Equivalence of Theorem 1 and the terminating case 
of Whipple's formula 



Invoke the terminating case of Whipple's summation formula (|B.9)1 with 

b = a + m + 3, c = a + 2 — m, (4.6) 

to evaluate the r.h.s. of (|4.5jl . and get 

~ , N Ar (-a - l) m ( a+2 - m - n )„ 
LJm,n) = 4" 



(a + l) m+2 (a + m + 3) n (a + 2 - m) n ' 



To identify this expression with (|1.6j) . decompose the ratio of the right sides of the two 
formulae into three groups and simplify using ()B.2|) - ()B.4|) : 



4"((a + d + 3)/2)„((a + rf + 4)/2)„ 1 



(a + l) m+ 2 (a + m + 3) n (a + d + 3) n a + V 
(-ar-l) m ((a + 2-m-n)/2) n 



(-a- l) d (a + 2-m) n ((d-a)/2) n 

Conversely, (jl.lljl implies ()B.9|) in the case when the parameters b and c can be presented 
in form ()4.6|) with m G Z + and Re a > — 1, that is if b — c is an even non-negative integer 
and Re (b + c) > 4. But both parts of (jB.9|) are rational functions of b and c . Hence the 
restrictions can be dropped. ■ 



4.3. Proof of Whipple's formula in the terminating case 

The l.h.s. of (|B.9|) is a rational function of a of the form 

. iF J ~n, n+1, a + 2 \ Q(a) 

\ a + m + 3, a + 2 — m J (a + m + 3) n (a — m + 2) n ' 

where is a polynomial of degree 2n (of course Q depends on m and n). We shall find 

all roots of Q(a) by using relation (|4.5jl and studying zeros and poles of L a (m,n). In this 
way we obtain a self-contained proof of Theorem 1 together with formula (jB.9|) . 



Lemma 4.3. If m—n >0 and even. Then (i) L a (m, n) = for a = 0, 2, . . . , m+n — 2 ; 
(ii) <5(m + n — 2s) = for s — 1, . . . , n. 

Proof, (i) It is enough to show that B a (m,n) = at the above values of a. In this 
case, \x — y\ a is a polynomial of degree < m + n. In its binomial expansion Ckix k y l , the 
exponents k, I satisfy at least one of the inequalities k < m, I < n . Therefore, all products 
(x k , p m {x)) (y l , p n (y)) are 0. 

(ii) Using the relation 

(-1 - a) m /(a + 2- m) n = (-1)" (-a - l) m _ n , 
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we get 

L a (m, n) = R(a) Q(a) a (a + 1) ... (a + m — n — 2), 

where R(a) is a rational function, which does not have poles and zeros at nonnegative 
integers. Now the claim follows from (i). □ 

In the next Lemma, we find n other roots of Q(a). 

Lemma 4.4. Let m > n. Then (i) L a (m, n) does not have poles at a = —m — n — 1 + 2s, 
s = 0, ...,n — 1; (ii) Q(— m — n — 1 + 2s) = for s = 0, . . . , n — 1. 

Proof, (ii) follows from (i), because the values a = —m — n — 1 + 2s, < s < n — 1, are 
uncompensated roots of the denominator in the r.h.s. in (J4.5|) . while L a (m,n) is regular at 
these values. 

To prove (i), consider analytical continuation of L a . Using the identity 

roc 

r a = (T(-a)y 1 e~ rt t- a - l dt (Re a < 0), 
Jo 

we obtain for — 1 < Re a < 

I oo 1 a; 

L a (m,n) = / r a_1 y e" te p m (x) y e ty p n (y) dy dx dt. 

^ ' 

Divide the t-domain into Jq 1 and J^ 00 ; the first part is a regular function of a for Re a < 0. 
Let us find asymptotics of the integrand as t — > oo in the second part and study its analytical 
continuation in Re a < — 1. 

The following formula of indefinite integration 

/ p(0e A5 d£ = e A «(A-V(0 - A-y (0 + A" V(0 - • • •) (4.8) 
holds for any polynomial So we have 



Applying ()4.8|) once again and using orthogonality, we get 

I e- tx p m (x) [ e ty Pn (y)dydx = £ £ pg>(0)rf?(0) + 0(e" 4 ) 

o j=0fc=0 

The coefficient at t _r is the residual of L a (m,n) at the pole a = — r, due to 

oo 1 

t~ a - l - r dt 



a + r 



By (jA.lOjl . these coefficients vanish for r — m + n + 1 — 2s, < s < n — 1. Hence there are 
no poles at the values stated in (i) of the Lemma. □ 
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The two Lemmas supply the complete list of 2n roots of the numerator Q(x). Thus we 
obtain: If m — n is even and > 0, then 

3 F 2 ( _n ' n+1 ' a + 2 ) C £S=I^ (^^T^)" (4 . 9) 
\a + m + 3, a + 2 — m) (a + m + 3) n (a — m + 2) n ' 

where C is independent of a. 

The left side is also a rational function of m with degrees 2n of the numerator and 
denominator, so the coefficient C is independent of m. Moreover, we can dismiss the above 
assumption of m — n being even and positive and allow m to be non-integral. To determine 
C (which now may depend only on n), let us take a — > oo. Then the l.h.s. of (14. 9j) is 
1 + o(l) , while the r.h.s. is C/4 n + o(l). Thus C = 4 n and (IB~9l) follows by substitution 
(HU). ■ 



4.4. Proof of Whipple's formula in the general case 

Fix the parameters b and c in Whipple's formula (|B.10J) . The difference of the l.h.s. and the 
r.h.s. of (jB.10|) is an entire function f(a). The motif of the proof is simple (cf. technique based 
on Carlson's theorem [JJ Chapter V]). The function f(a) has sufficiently many known zeros 
and a uniformly controlled growth as \a\ — > oo. From this we will conclude that f(a) = 0. 

At this time, the proof is not self-contained, though marginally. I have to recourse to 
additional facts from the theory of generalized hypergeometric series: Karlsson's formula 
and Dixon's formula. 

Lemma 4.5. (i) // Re (b + c) > 1 then the right-hand side of W.10\) has estimation 

0( yr]Ima| ^ + | a |)l-R*(6+c)) ag | a | _> qq. 

(ii) // b, c are real and b + c> 1 then the left-hand side of WJ(\) is 0{e^ a \ (1 + |a|) 1_ ( 6+c )) 
as \a\ — > oo. 

Proof, (i) follows from the inequality 

(T(z)T(r - z))- 1 < C , (r)e" |Im2| (l + |z|) 1 - Rer , (4.10) 

which is a consequence of the formula r(;?)r(l — z) — ir/ sin(7rz) and the estimation 
|r(x)/r(x + r)\ = 0(x~ Rer ), x — > oo, I argx| < n. 

(ii) The series has the form 

~ (ojn (1 - a) n 

t/ n n\((b + c+l)/2) n > 

where 

((6 + c-l)/2) B ((6 + c+l)/2) n 
/in = 7T\ / x = 0(1) as n > oo. 

[b)n {C)n 

Next, (a) n (1 — a) n is a series with positive coefficients in t = (a—1/2) 2 . So maximum of our 
3 F 2 on the circle \a — 1/2 1 = R is bounded, up to a constant, by the respective maximum of 



l F 1 



a. 



(b + c+l)/2j ' r(^±fti - a) r(^±|=i + a) ' 
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Finally, by lOTO) . the r.h.s. is 0(e^ Ima l (1 + lal) 1 '^). □ 

If in (ii) we had e^' 1 ™ ' instead of e w ^ , then the proof could be easily completed via 
the maximum principle applied to the function f(a) cosecna, which is entire. (No poles on 
the terminating case is already established). 

In reality, we have only found that In \ f{a)\ < (it + o(l))|a|. Consider the asymptotical 
density of zeros of /, 



v= lim N(R)/R, N(R) = #{a: f(a) = and \a\ < R } 

R—>oo 

It follows from Jensen's formula that the estimation v > 2 would imply / = 0. Thus far 
we only know that f\Z = 0, which gives v > 2 but leaves a possibility of the equality. To 
complete the proof, we shall present another infinite sequence of zeros. Notice that the r.h.s. 
of (jB.10|) vanishes when a = — b — 2k, k G 1L. Show that those are also zeros of the l.h.s. 
(By symmetry between b and c, yet more zeros will be found.) 

Lemma 4.6. Let n be a nonnegative even integer, and b, c arbitrary, except nonnegative 
integers, and denote d — (b + c — l)/2. Then 



n, b + n + 1, cf 



b, c 



0. (4.11) 



Proof. Through Karlsson's formula (jB.12|) and Gauss' formula ()B.7|) we transform (|4.11|) 
into the terminating 3 F 2 with parameters (— n — 1, —b — n, d\b, —d — n), times a nonsingu- 
lar factor. Such an 3 F 2 falls under Dixon's formula [SI 3.1], [UJ 4.4 (5)], the sum being a ratio 
of products of T functions. The only singular factor, T(—n), appears in the denominator. 
(The factor T((l — n)/2) is safe, because n is even.) Therefore, the ratio equals to 0. ■ 



§ 5. Proof of Theorem 4 

The double integral L'(m,n) will be expressed via a double sum using the Maclaurin ex- 
pansions of the shifted Legendre polynomials. We rearrange the double sum and show that 
new partial single sums satisfy certain recurrence relation, from which the theorem follows. 

We set aside the simple case m = n = and prove formulae ()1.15j) and ()1.16|) assuming 
m = max(m, n) > 0. 

5.1. Reduction of the integral B' to a finite sum 

Below any function means a function as good as a polynomial. 

Denote the average value of a function / on an interval (a, b) by 

</>£=(& -a) -1 [ b f(t)dt. 
Introduce two operators, M and N, as follows: 

Mf(x) = {f)l (5.1) 
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Nf(x)= / (fr t dt. (5.2) 



o 



Clearly, if / is a polynomial of degree n then Mf is a polynomial of degree n , and Nf a 
polynomial of degree n+1. 

Lemma 5.1. Let f , g be any functions, and Gix) = \ giy)dy. Then 

Jo 

1 X 1 1 

J S{x) J\n(x ~y)g(y)dydx = - J M[f ■ G}(x) dx - J f(x) ■ Ng(x)dx. (5.3) 



Proof. Integration by parts w.r.t. y yields 

rx 

/ \n(x — y) g(y) dy = \nxG(x) — Ng(x). 
Jo 

The assertion of Lemma follows by second integration 

f 1 lux f{x) G(x) dx = - f 1 M[f ■ G)(x) dx, 
Jo Jo 

where we remember that \im x ^oG(x) Inx = 0. □ 

Definition of the sums S(m,n,r). Let < r < (m + n). The quantities S(m,n,r) are 
the coefficients of the generating function 



m+n 



P m (x)q n+1 (x) = S{m,n,r)x r+1 , (5.4) 

r=0 

where q n+ \ are the integrated shifted Legendre polynomials — see (|A.12|) . The explicit 
formula is (cf. (TOTl) and (TQ2l) ) 

S( m ,n,r) = (-1)-" if j^^ -, (",+,*»' (5.5) 

j= 7^ in (m-j)!j! 2 (n-k)\k\{k + l)\ 

where k = r — j, jmin = max(0, r — n), jmax = min(m, r). 
The promised double sum and its rearrangement are as follows 



j=0 k=0 



(m-j)\j\ 2 (n-k)\k\(k + l)\ (j + k + 2) 



m+n S{m, n, r 



(5.6) 



V"") ' ) def T r, \ 

w = K{m > n) - 



Lemma 5.2. For any m, n 

Pm{x)p n {y) \n{x-y)dydx = -K(m,n) - (p m , Np n ). (5.7) 



0<y<x<l 
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Proof. Evaluate the functional u — > (Mu, 1) , which appears in (|5.3jl . in terms of the 
Maclaurin expansion of the function u(x) 

The operator M divides the r-th coefficient by r + 1. Thus 

(Mu, l> = J> r (r + 1)" 2 . 

Now take in (|5.3|) / = p m , g = p n , hence u = p m q n +i ■ Replacing u r+ i by S(m,n,r), we 
arrive at (|5.7j) . □ 

The scalar product (p m , Np n ) can be found without difficulty, if m > n. 

Lemma 5.3. (i) If m > n + 2 , then (p m , Np-n) = 0. 
(ii) Nonvanishing scalar products with m > n are 

( Pn , Np n ) = [(2n + l)(2n + 2)]-\ (5.9) 

^ NPn) = 2(2n+l) + (2n + 3y ^ 
where H r denotes the r-th harmonic sum 

H r = 1 + 2- 1 + ...+ r _1 . (5.11) 

Proof, (i) follows from orthogonality, since deg Np n = n + 1 < in. 

(ii) Apply the operator iV to a Maclaurin series (|5.8|) . If U = J u, then 

U{ X )-U{t) = ^_ r 
r — t ^ r + 1 



Therefore, 



Nu(x) = J2u r ^\x r+1 . (5.12) 

r>0 r + 1 



Considering the expansion of Np n , we obtain (in notation (jA.lljO 

iVp n (x) = < j -I + +On-l— X 

n + 1 n 
„(n) ^ / „(«) J, „(f»+l) /, 



^( n+1 )n + l I />) n n (n+1) n + l 

a n+ i '* ^ 1 \ On " a n+1 "ti 

Since n -1 a^l x /affi = — 1/2 is independent of n, the coefficient at p n in iVp n equals 
(-1/2) (h n - h n+ i) = (2n + 2)" 1 . Recalling the value (IA.8I) of (p n ,Pn), we obtain 
The coefficient at p n+ i in Afy>„ and consequently the value (|5.10|) also follow by ()A.11|) and 
(TOl). □ 
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5.2. Recurrence for S(m,n,r). 
Lemma 5.4. The following expression is symmetric in m and n (m,n > 1) 

D(m,n,r) = (n + 1) S(m, n, r) — (n — 1) S(m — 1, n — 1, r) (5.13) 

Proof. The generating function D mn (x) for the quantities D(m,n,-) is 

m+n 

D(m,n,r) x r+1 = (n + 1) p m (x) q n+1 (x) - (n - l)p m _i(x) q n {x). 

r=0 

By (EH, 

n+1 n-1 

u mn — ~ — Vm {Pn+1 ~ Pn~l ~ ~ 7 Pm-1 \Pn ~ Pn-2) ■ 

In + 1 In — 1 

Apply (JA.9j) to Pn+i and p n and find 

n / \ _ ~ n — 1 _ n-1 

Dmn\ x ) — x PmPn ~ PmPn-1 ~ % Pm-lPn-1 + Pm-lPn-2- 

n n 



The first term is explicitly symmetric. Apply (jA.9|) to p m in the 2nd term and obtain a 
symmetric formula for D mn (x) — xp m p n : 

f n — 1 m — 1\ f m — 1 n — 1 \ 

1 H 1 Xp m _ip n _i + Pm-2Pn-l H Pm-lPn-2 • 

Vram/ V m n J 

□ 

Corollary. K(m, n) satisfy the recurrence relation 
(n + 1) K(m, n) — (m + 1) i^(n, m) 

= (n — 1) _fT(m — 1, n — 1) — (m — 1) K(n — l,m — 1). 



(5.14) 



5.3. Completion of proof 

Lemma 5.5. (i) If \m — n\ > 2 then K(m,ri) + K(n,m) = 0. 
(ii) In £/ie cases when \m — n\, we have 

2K{n,n) = [n(2n + l)(2n + 2)] _1 . (5.15) 

if (n + 1, n) + if (n, n + 1) = -[(2n + l)(2n + 2)(2n + 3)]" 1 . (5.16) 

Proof. From proof of Lemma 15.2.1 we see that 

K{m,n) + K(n,m) = (M[p m q n+l + q m+1 p n ], 1) = / gm+1 ^ 9n+1 X dx. 

o 

If n < m — 1 then g m+ i is orthogonal to any polynomial of degree n (see (|A.13)l ). in 
particular to x _1 g n+ i(x). This proves (i). 
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(ii) Set 

w n +i(x) = I x~ l q n+ i(x)dx. 



o 



Integrating by parts and remembering that g n+ i(0) = q n+ i(l) = for n > 0, we get 



K(m,n) + K(n,m) — — I p m (x) w n+ i(x) dx. 



o 



Find the (n + l)-th and n-th Legendre coefficients of u> n +i by considering the respective 
Maclaurin coefficients. The calculation is very similar to that in proof of Lemma 15.3.1 



l _n x n+l _|_ a n-l 

(n + l) 2 n 



W n+ l{X) = — ; 1N9 X + — —X 



Pn + i(x) ( a { :l, 1 aj w+1) , 

(r,\ _9 (Vi-HI , 1\5 PnW 



a (n+l) (n+l) 2 ^ fl (n) n 2 Q (n+1) (n + l) 

The coefficient at p n equals (—1/2) [n^ 1 — (n + 1) _1 ] = — [2n(n + The coefficient at 

p ra+ i is [(2n + l)(2n + 2)] _1 . The proof is completed by using (|A.8|) . as in Lemma 15.3.1 □ 

Lemma 5.6. If m > 2, then 

K(m } 0) = (-l) m+1 [m(m + 2)(m 2 - l)]" 1 . (5.17) 

Proof. By (fA~9|) . 

2p m (x)gi(x) = (x + l)p m (x) = ^ + ^ Pm+l(g) + Pm(&) + ~ : p m -l(g). 

2m + 1 2m + 1 

Applying the operator M and integrating, we get 

2K(m 0) = j /m+j_ gm+2 (x) + U^) + ^ 

7 V 2m + 1 x x 2m + 1 x I 

o v 7 

Now use (lA~T4l) . and (l5~T7f) follows . □ 

Final evaluation of L'(m,n) requires separate treatment of the cases m — n > 2, m = n, 
and m — n — 1. 

Case I : m — n > 2. In this case, K(m,n) = —K(n,m) by Lemma 15.5.1 (i), and 1)5.14)1 

becomes 

(m + n + 2) K(m, n) = (m + n — 2) f^(m — 1, n — 1). 
The value K(m, 0) is known from (J5.17)) . By induction we find 

K(m, n) = — [(m + n)(m + n + 2)((m — n) 2 )] -1 . 

Since the second term in the r.h.s. of ()5.7)) vanishes — see Lemma f5.3.l (i) — we obtain the 
answer (jl.15)) . 
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Case II : m = n. Using ()5.15|) and ()5.9|) . we find the value of ()5.7|) : 

1 1 1 



2n (2n + 1) (2n + 2) (2n + 1) (2n + 2) 2n(2n + 2)(-l)' 
which agrees with (jl.l5|) . 

Case HI : m — n + 1. Denote if(n + l,n) = A(n). The recurrence (|5.14|) with the 
nonhomogeneous symmetry relation f)5.16j) leads to the nonhomogeneous linear difference 
equation 

2n + 5 

(2n + 3 )A„-(2„-l)A„- 1= (2?t _ 1)(2n+1)3 . 

A solution of the homogeneous part is = [(2n + l)(2n + 3)] _1 . Set \ n = fi n . Then 
fi n satisfies the equation 

2n + 5 1 / 1 1 

Hn — — J^. TTTT; i \ — o 7. 7 + 



(2n- l)(2n + 2) 2 2 V2n - 1 2n + 3 n + 1 
By telescoping, 



Hn = to - * + #2°n+l - n + - ((2n + - (2n + l)' 1 ) . (5.18) 



6 + 2n+1 ~ 2 ^ n+1 + 2 

Explicit calculation shows that /x = 3if(l, 0) = —1/12. Substituting (|5.18|) together with 
(EHUl) to (|577jl . we obtain (fTTHl). 

The proof is complete. ■ 
Excercise. Let m, n, / be nonnegative integers. Evaluate the sum: 

T{m,n,l) = ^ . 



Appendix A. Legendre polynomials P n (x) 

Definition. P n {x) is the polynomial of degree n such that: (i) -P n (l) = 1, and (ii) for 
any polynomial f{x) of degree < n 

i 

(/, P n ) = J f(x)P n (x)dx = 0. 
-l 

Rodrigues formula 

1 d n 

p( x ) = — (x 2 - 1)". (A.l) 

Symmetry 

P n (-x) = (-l) n P n (x). (A.2) 
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Li2 norm 

(P n , P n ) = 2{2n + iy l . (A.3) 

Recurrence relations 

(n + l)P n+1 (x) = (2n + l)xP n (x) - nP n _i(z), (A.4) 

P' n+l {x) = (2n + 1) P„(x) + ^(x). (A.5) 
Shifted Legendre polynomials [1, 22.2.11] 

p n (x) = P n (x), x = 2x - 1, (A. 6) 

form an orthogonal basis on [0, 1]. When working with p n 's, we denote by (•, •) the integral 
scalar product on [0, 1]. For reference, here are the shifted versions of the above formulae: 

p n ( x ) = (d)- 1 (d/dx) n ( x n {x - l) n ) , (A.7) 

(p n , p n ) = {2n + I)' 1 , (A.8) 

(n + l)p n+1 (x) = (2n + l)xp n (x) - np n -i(x). (A. 9) 
A bilinear identity for derivatives. If r = \m — n\ + 2s, s > 0, then 

E (-l) j ^ ) (0)rf' ) (0) = 0. (A.10) 

j+k=r+l 

Proof. The same combination of derivatives taken at 1 has opposite sign. So, by integration 
by parts, the l.h.s. equals to 

(±1/2) ( <p m , p^) - {pt +l \p n )\ 

which is 0, because, e.g., deg(p^ +1 )) = max(0, n — 1 — 2s) < n. □ 
Maclaurin expansion 

f ^ V- (n) k in) (-l)"~ fc (n + fc)! 

Pn{X) = E °fc X > = —7 aim • ( A - U ) 



Integrated shifted Legendre polynomials g n+ i(x) 

,(n) 
+ 



q n+ i(x) = / p^dt = '£-±-x k+1 . (A.12) 
Jo fc=0 fc + 1 



These polynomials are used in Sect. 5. The notation is local to this work. 
Decomposition (Cf. (|A.5)1 .) For n > 0, 

2 (2n + 1) g„+i = p n+ i - p n -i- (A. 13) 

Integral (g n+ i, x^ 1 ). If n > 0, then 





y x y ^ v 7 n n + 1 v y 



n(n + 1) 
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Proof. Set w(x) = x(x — 1). By the Rodrigues formula q n+ i = (n\) 1 d n 1 w n . Integrating 
by parts n times, we get 

[n — \)\ [ w n 

■ ax. 



J x n n\ x3 +1 n\ J 



X'' 



1 



When e — > + , the only survivor is n 1 J(x — l) n dx = (— l) n [n(n + 1)] □ 



o 



Appendix B. Hypergeometric series and identities 

Rising factorial (Pohhammer's symbol). Let n G Z, a £ C. By definition, (a) n = 
r(a + n)/r(a). In particular, (l) n = n\. Alternative definition: 

(i) (o)o = i ; 

(ii) (a) n = a(a + l) ■ ... • (a + n- 1) if n > 0; (B.l) 

(iii) (a)„ = [(a + n)_ n ] _1 if n < 0. 

Three useful formulae: 

{a) m+n = (a) m (a + m)„, (B.2) 

(a) n = (-l) n (-a-n + l) n , (B.3) 

(a) n (a + l/2) n = 4-" (2a) 2n , (B.4) 
Hypergeometric series. The series depending on p + g + 1 parameters 

/V . . . , Op \ = ~ (ai) r , . . (a p ) n „ (B 5) 

is called the (generalized) hypergeometric series. The ordinary (Gauss') series is 2 -Pi- If one 
of the a,j 's in f)B.5|) is a negative integer — r, then the series turns to a polynomial of degree 
r in 2 and is called terminating. In the frequent case z = 1, the argument z of p F g is 
omitted. 

Euler's integral for 2 F X [1, (15.3.1)], [6., 2.1.3 (10)], 5, 1.5]. 
If Rec > Re 6 > and z ^ [1, 00), then 

a Fi ^ ;z) = J t b ~ l (1 - tf- b - 1 (1 - tz)- a dt. (B.6) 
^ ' 

Gauss' formula for 2 F 1 (. . . ; 1) (15.1.20)], 2.1.3 (14)], 5, 1.3] 

^ fa, b\ T(c)T(c — a — b) , . . . 

»M c) = rc-.)r(c-6 — »)>»)■ 
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Thomae's transformation formula for 3 F 2 O (2)], 3.2], jHJ (3.1.2)] If s = 
e + f — a — b — c, then 

fa, b, c\ = T(e)T(f)T(s) fe-a, f -a, s\ 

3 2 \ e, f J T(a)T(s + b)T(s + c) 3 2 ^ s + b, s + c )' 1 " ' 

Whipple's summation formula for 3 F 2 [14, (34)], 3.4], 4.4 (7)] . Terminating 
case: for a nonnegative integer n, 

3F J-n, » + l (H C -l)M 4 »P|M 

V b, C J {0)n{c)n 

General case: let a e C, a' = 1 — a, and Re (6 + c) > 1. Then 

„ fa, a', (6 + c-l)/2\ 4tt -2-( 6+c )r(6)r(c) 



3-^2 



ft, c ; r(^)r(^)r(^)r(^ 



Dougall's formula for 5 F 4 4.4], [UJ 4.5 (6)] . Suppose that 



1 + a = b + b' = c + c = d + d! = e + e' 



and b — b' = 1. Then 

'a, b, c, d, e\ r(c') T{d') r(e') T(l + a - c - d - e) 



5^4 



b', d, d', e'J T(l + a) T(c' - d) T(d' - e) T(e' - c) 
Particular 3 F 2 case of Karlsson's reduction formula [9 



3 r 2 



—b — n, b + n+1, d\ 

b, c ) 

V*k fn + 1\ {-b - n) k {d) k f-b-n + k, d + k^ 
= u m" f „\. 2i<1 



k=0 



k ) (b) k (c) k \ c + k 



(B.10) 



(B.ll) 



(B.12) 
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